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Abstract: 

In the latter half of the last century, it became evident that there exists an ever increasing number 
of different states of the so-called elementary particles. The usual reductionist approach to 
this problem was to search for a simpler infrastructure, culminating in the formulation of the 
quark model and quantum chromodynamics. In a complementary, completely novel approach, 
Hagedorn suggested that the mass distribution of the produced particles follows a self-similar 
composition pattern, predicting an unbounded number of states of increasing mass. He then 
concluded that such a growth would lead to a limiting temperature for strongly interacting 
matter. We discuss the conceptual basis for this approach, its relation to critical behavior, and 
its subsequent applications in different areas of high energy physics. 


* To appear in Melting Hadrons, Boiling Quarks, R. Hagedorn and J. Rafelski (Ed.), Springer 
Verlag 2015, open access. 



A prophet is not without honour, 
but in his own country. 

The New Testament, Mark 6,4. 


1 Rolf Hagedorn 

The development of physics is the achievement of physicists, of humans, persisting against 
often considerable odds. Even in physics, fashion rather than fact is frequently what determines 
the judgement and recognition. 

When Rolf Hagedorn carried out his main work, now quite generally recognized as truly pio¬ 
neering, much of the theoretical community not only ignored it, but even considered it to be 
nonsense. “Hagedorn ist ein Narr”, he is a fool, was a summary of many leading German theo¬ 
rists of his time. When in the 1990’s the question was brought up whether he could be proposed 
for the Max Planck Medal, the highest honor of the German physics community, even then, 
when his achievements were already known world-wide, the answer was still “proposed, yes...” 
At the time Hagedorn carried out his seminal research, much of theoretical physics was ideo¬ 
logically fixed on “causality, unitarity, Poincare invariance”: from these three concepts, from 
axiomatic quantum field theory, all that is relevant to physics must arise. Those who thought 
that science should progress instead by comparison to experiment were derogated as “fitters and 
plotters”. Galileo was almost forgotten... Nevertheless, one of the great Austrian theorists of 
the time, Walter Thirring, himself probably closer to the fundamentalists, noted: “If you want 
to do something really new, you first have to have a new idea”. Hagedorn did. 

He had a number of odds to overcome. He had studied physics in Gottingen under Richard 
Becker, where he developed a life-long love for thermodynamics. When he took a position at 
CERN, shortly after completing his doctorate, it was to perform calulations for the planning and 
construction of the proton synchrotron. When that was finished, he shifted to the study of mul¬ 
tihadron production in proton-proton collisions and to modelling the results of these reactions. 
It took a while before various members of the community, including some of the CERN The¬ 
ory Division, were willing to accept the significance of his work. This was not made easier by 
Hagedorn’s strongly focussed region of interest, but eventually it became generally recognized 
that here was someone who, in this perhaps a little similar to John Bell, was developing truly 
novel ideas which at first sight seemed quite specific, but which eventually turned out to have a 
lasting impact also on physics well outside its regions of origin. 

We find that Rolf Hagedorn’s work centers on two themes: 

• the statistical bootstrap model, a self-similar scheme for the composition and decay of 
hadrons and their resonances; for Hagedorn, these were the “fireballs”. 

• the application of the resulting resonance spectrum in an ideal gas containing all possible 
hadrons and hadron resonances, and to the construction of hadron production models 
based on such a thermal input. 

We will address these topics in the first two sections, and then turn to their role both in the 
thermodynamics of strongly interacting matter and in the description of hadron production in 
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elementary as well as nuclear collisions. Our aim here is to provide a general overview of 
Hagedom’s scientific achievements; other aspects will be covered in other chapters of this book. 
Some of what we will say transcends Hagedom’s life. But then, to paraphrase Shakespeare, we 
have come to praise Hagedom, not to bury him; we want to show that his ideas are still important 
and very much alive. 


2 The Statistical Bootstrap 

Around 1950, the world still seemed in order for those looking for the ultimate constituents 
of matter in the universe. Dalton’s atoms had been found to be not really atomos, indivisible; 
Rutherford’s model of the atom had made them little planetary systems, with the nucleus as the 
sun and the electrons as encircling planets. The nuclei in turn consisted of positively charged 
protons and neutral neutrons as the essential mass carriers. With an equal number of protons and 
electrons, the resulting atoms were electrically neutral, and the states obtained by considering 
the different possible nucleus compositions reproduced the periodic table of elements. So for a 
short time, the Greek dream of obtaining the entire complex world by combining three simple 
elementary particles in different ways seemed finally feasible: protons, neutrons and electrons 
were the building blocks of our universe. 

But there were those who rediscovered an old problem, first formulated by the Roman philoso¬ 
pher Lucretius: if your elementary particles, in our case the protons and neutrons, have a size 
and a mass, as both evidently did, it was natural to ask what they are made of. An obvious 
way to find out is to hit them against each other and look at the pieces. And it turned out that 
there were lots of fragments, the more the harder the collision. But they were not really pieces, 
since the debris found after a proton-proton collision still also contained the two initial protons. 
Moreover, the additional fragments, mesons and baryons, were in almost all ways as elemen¬ 
tary as protons and neutrons. The study of such collisions was taken up by more and more 
laboratories and at ever higher collision energies. As a consequence, the number of different 
“elementary” particles grew by leaps and bounds, from tens to twenties to hundreds. The latest 
compilation of the Particle Data Group contains over a thousand. 

What to do? One approach was in principle obvious: just as Dalton’s atoms could be con¬ 
structed from simpler, more elementary constituents, so one had to find a way of reproducing 
all the hadrons, the particles formed in strong interaction collisions, in terms of fewer and more 
elementary building blocks. This conceptually straight-forward problem was, however, far from 
easy: a simply additive composition was not possible. Nevertheless, in the late 1960’s the quark 
model appeared, in which three quarks and their antiquarks were found to produce in a non- 
Abelian composition all the observed states and more, predicted and found. Not much later, 
quantum chromodynamics (QCD) appeared as the quantum field theory governing strong in¬ 
teractions; moreover, it kept the quarks inherently confined, without individual existence: they 
occurred only as quark-antiquark pair (a meson) or as a three-quark state (a baryon). And 
wealth of subsequent experiments confirmed QCD as the basic theory of strong interactions. 
The conventional reductionist approach had triumphed once more. 

Let us, however, return to the time when physics was confronted by all those elementary par¬ 
ticles, challenging its practioners to find a way out. At this point, in the mid 1960’s, Rolf 
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Hagedom came up with a truly novel idea lUl - S . He was not so mueh worried about the spe- 
eifie properties of the partieles. He just imagined that a heavy partiele was somehow eomposed 
out of lighter ones, and these again in turn of still lighter ones, and so on, until one reaehed the 
pion as the lightest hadron. And by eombining heavy ones, you would get still heavier ones, 
again: and so on. The crucial input was that the eomposition law should be the same at each 
stage. Today we call that self-similarity, and it had been around in various forms for many years. 
A partieularly elegant formulation was written a hundred years before Hagedorn by the English 
mathematieian Augustus de Morgan, the first president of the London Mathematieal Soeiety: 

Great fieas have little heas upon their backs to bite’em, 
and little fleas have lesser still, and so ad infinitum. 

And the great fleas themselves, in turn, have greater fleas to go on, 
while these again have greater still, and greater still, and so on. 

Hagedom proposed that “a fireball eonsists of fireballs, whieh in turn eonsist of fireballs, and 
so on...” The eoneept later reappeared in various forms in geometry; in 1915, it led to the 
eelebrated triangle devised by the Polish mathematieian Waelaw Sierpinski: “ a triangle eonsists 
of triangles, whieh in turn eonsist of triangles, and so on...”, in the words of Hagedom. Still 
later, shortly after Hagedom’s proposal, the French mathematician Benoit Mandelbroit initiated 
the study of such fractal behaviour as a new field of mathematies. 



Hagedorn had reealled a similar problem in number theory: how many ways are there of de- 
eomposing an integer into integers? This was something already adressed in 1753 by Leonhard 
Euler, and more than a eentury later by the mathematieian E. Sehroder in Germany. Finally G. 
H. Hardy and S. Ramanujan in England provided an asymptotie solution. Let us here, however, 
eonsider a simplified, easily solvable version of the problem in whieh we eount all possible 
different ordered arrangements p{n) of an integer n. So we have 


1=1 

p{l) = 1 =2"-i 

2=2,1-1-1 

p(2) =2 = 2"-i 

3=3, 2-1-1, 1-1-2, l-i-l-i-l 

p(3) = 4 = 2"-! 

4= 4, 3-1-1, 1-1-3, 2-1-2, 2-I-1-I-1, l-i-2-i-l, 1-I-1-I-2, l-i-l-i-l-i-l 

p{A) = S = 2"-'^ 

and so on. In other words, there are 


p(n) = 2"'i = - 

(1) 


4 







ways of partitioning an integer n into ordered partitions: p{n) grows exponentially in n. In 
this particular case, the solution could be found simply by induction. But there is another 
way of getting it, more in line with Hagedom’s thinking: “large integers consist of smaller 
integers, which in turn consist of still smaller integers, and so on...” This can be formulated as 
an equation, 

p(n) = 5(n-l)+ ^ n (2) 

It is quite evident here that the form of the partition number p (n) is determined by a convolution 
of many similar partitions of smaller n. The solution of the equation is in fact just the number 
of partitions of n that we had obtained above. 


p{n)=zp{n) (3) 

up to a normalization constant of order unity (for the present case, it turns out that z ~ 1.25). 
For Hagedorn, eq. Q expressed the idea that the structure of p(«) was determined by the 
structure of p(n) - we now call this self-similar. He instead thought of the legendary Baron von 
Miinchhausen, who had extracted himself from a swamp by pulling on his own bootstraps. So 
for him, eq. Q became a bootstrap equation. 

The problem Hagedorn had in mind was, of course, considerably more complex. His heavy 
resonance was not simply a sum of lighter ones at rest, but it was a system of lighter resonances 
in motion, with the requirement that the total energy of this system added up to the mass of the 
heavy one. And similarly, the masses of the lighter ones were the result of still ligher ones in 
motion. The bootstrap equation for such a situation becomes 


p{m,Vo) = d{m-mo) 

N 


1 

m 


Vo 



{nii) (Ppi] d^{Y.iPi-p), 


(4) 


where the first term corresponds to the case of just one lightest possible particle, a “pion”. The 
factor Vq, the so-called composition volume, specified the size of the overall system, an intrinsic 
fireball size. Since the mass of any resonance in the composition chain is thus determined by 
the sum over phase spaces containing lighter ones, whose mass is specified in the same way, 
Hagedorn called this form of bootstrap “statistical”. 

After a number of numerical attempts by others, W. Nahm [|7]| solved the statistical bootstrap 
equation analytically, obtaining 


p(m,Vb) = const, m ^exp{m/r/f}. 


(5) 


So even though the partitioning now was not just additive in masses, but included the kinetic 
energy of the moving constituents, the increase was again exponential in mass. The coefficient 
of the increase, 7)^ \ is determined by the equation 

^^{molTHfK 2 imolTH) = 21n2 - 1, (6) 

in terms of two parameters Vq and mo. Hagedorn assumed that the composition volume Vq, 
specifying the intrinsic range of strong interactions, was determined by the inverse pion mass 
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as scale, Vq — (4^/3)m;j^. This leads to a scale factor ~ 150 MeV. It should be emphasized, 
however, that this is just one possible way to proceed. In the limit mo —)■ 0, eq. ® gives 


r^ = [;r2(21n2-l)]i/3yo ~ 1/r,, 


(7) 


where Vq = (47r/3)r| and denotes the range of strong interactions. With r/j ~ 1 fm, we thus 
have Th — 200 MeV. From this it is evident that the exponential increase persists also in the 
chiral limit m-ji —0 and is in fact only weakly dependent on mo, provided the strong interaction 
scale Vq is kept fixed. 

The weights p(m) determine the composition as well as the decay of “resonances”, of fireballs. 
The basis of the entire formalism, the self-similarity postulate - here in the form of the statis¬ 
tical bootstrap condition - results in an unending sequence of ever-heavier fireballs and in an 
exponentially growing number of different states of a given mass m. 

Before we turn to the implications of such a pattern in thermodynamics, we note that not long 
after Hagedom’s seminal paper, it was found that a rather different approach, the dual resonance 
model [ISUTOll led to very much the same exponential increase in the number of states. In this 
model, any scattering amplitude, from an initial two to a final n hadrons, was assumed to be 
determined by the resonance poles in the different kinematic channels. This resulted structurally 
again in a partition problem of the same type, and again the solution was that the number 
of possible resonance states of mass m must grow exponentially in m, with an inverse scale 
factor of the same size as obtained above, some 200 MeV. Needless to say, this unexpected 
support from the forefront of theoretical hadron dynamics considerably enhanced the interest in 
Hagedom’s work. 


3 The Limiting Temperature of Hadronic Matter 


Consider a relativistic ideal gas of identical neutral scalar particles of mass mq contained in a 
box of volume V, assuming Boltzmann statistics. The grand canonical partition function of this 
system is given by 



N 


( 8 ) 


leading to 


\nnT,V) = -^K2{'^). 


VTml 


(9) 




For temperatures T mq, the energy density of the system becomes 



( 10 ) 


the particle density 
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( 11 ) 



and so the average energy per partiele is given by 


(0 ~ 3 r. 


( 12 ) 


The important feature to learn from these relations is that, in the ease of an ideal gas of one 
speeies of elementary partieles, an inerease of the energy of the system has three eonsequenees: 
it leads to 

• a higher temperature, 

• more eonstituents, and 

• more energetie eonstituents. 

If we now consider an interacting gas of such basic hadrons and postulate that the essential form 
of the interaction is resonance formation, then we can approximate the interacting medium as 
a non-interacting gas of all possible resonance species ifTTlfT^ . The partition function of this 
resonance gas is 



where the sum begins with the stable ground state mo and then includes the possible resonances 


niij = 1,2,... with weights p{mi) relative to mo. Clearly the crucial question here is how 
to specify p{mi), how many states there are of mass m,. It is only at this point that hadron 
dynamics enters, and it is here that Hagedom introduced the result obtained in his statistical 
bootstrap model. 

As we had seen above in eq. dS]), the density of states then increases exponentially in m, with 
a coefficient ^ determined by eq. ® in terms of two parameters Vo and niQ. If we replace 
the sum in the resonance gas partition function (fT3l) by an integral and insert the exponentially 
growing mass spectrum ® , eq. (fTSl) becomes 




(14) 


Evidently, the result is divergent for all T > Th'. in other words, Th is the highest posssible 
temperature of hadronic matter. Moreover, if we compare such a system with the ideal gas of 
only basic particles ( a “pion” gas), we find 


pion gas 


resonance gas 



tires ~ £ 

(Ores ~ const. 


Here n denotes the average number density of constituents, (O the average energy of a con¬ 
stituent. In contrast to to the pion gas, an increase of energy now leads to 
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• a fixed temperature limit, T ^Th, 

• the momenta of the constituents do not continue to increase, and 

• more and more species of ever heavier particles appear. 

We thus obtain a new, non-kinetic way to use energy, increasing the number of species and 
their masses, not the momentum per particle. Temperature is a measure of the momentum of 
the constituents, and if that cannot continue to increase, there is a highest possible, a “limiting” 
temperature for hadronic systems. 

Hagedom originally interpreted Th as the ultimate temperature of strongly interacting matter. 
It is clear today that Th signals the transition from hadronic matter to a quark-gluon plasma. 
Hadron physics alone can only specify its inherent limit; to go beyond this limit, we need more 
information: we need QCD. 

As seen in eq. ([5]), the solution of the statistical bootstrap equation has the general form 

p(m, Vb) ~ "i~^exp{m/r//}, (15) 

with some constant a; the exact solution of eq. dH) by Nahm gave a = 3. It is possible, however, 
to consider variations of the bootstrap model which lead to different a, but always retain the 
exponential increase in m. While the exponential form makes Th the upper limit of permissible 
temperatures, the power law coefficient a determines the behavior of the system at T = Th- 
For a = 3, the partition function (fT4l) itself exists at that point, while the energy density as 
first derivative in temperature diverges there. This is what made Hagedom conclude that Th is 
indeed the highest possible temperature of matter: it would require an infinite energy to reach 
it. 

Only a few years later it was, however, pointed out by N. Cabibbo and G. Parisi ifTSl that larger 
a shifted the divergence at T = Th to ever higher derivatives. In particular, for 4 > a > 3, the 
energy density would remain finite at that point, shifting the divergence to the specific heat as 
next higher derivative. Such critical behavior was in fact quite conventional in thermodynamics: 
it signalled a phase transition leading to the onset of a new state of matter. By that time, the 
quark model and quantum chromodynamics as fundamental theory of strong interactions had 
appeared and suggested the existence of a quark-gluon plasma as the relevant state of matter at 
extreme temperature or density. It was therefore natural to interpret the Hagedom temperature 
Th as the critical transition temperature from hadronic matter to such a plasma. This interpre¬ 
tation is moreover corroborated by a calculation of the critical exponents IfTSlI governing the 
singular behavior of the resonance gas thermodynamics based on a spectrum of the form (fT5l) . 

It should be noted, however, that in some sense Th did remain the highest possible temperature 
of matter as we know it. Our matter exists in the physical vacuum and is constmcted out of 
fundamental building blocks which in turn have an independent existence in this vacuum. Our 
matter ultimately consists of and can be broken up into nucleons; we can isolate and study a 
single nucleon. The quark-gluon plasma, on the other hand, has its own ground state, distinct 
from the physical vacuum, and its constituents can exist only in a dense medium of other quarks 
- we can never isolate and study a single quark. 

That does not mean, however, that quarks are eternally confined to a given part of space. Let 
us start with atomic matter and compress that to form nuclear matter, as it exists in heavy nu¬ 
clei. At this stage, we have nucleons existing in the physical vacuum. Each nucleon consists of 
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three quarks, and they are confined to remain close to each other; there is no way to break up a 
given nucleon into its quark constituents. But if we now continue to compress, then eventually 
the nucleons will then penetrate each other, until we reach a dense medium of quarks. Now 
each quark finds in its immediate neighborhood many other quarks besides those which were 
with it in the nucleon stage. It is therefore no longer possible to partition quarks into nucleons; 
the medium consists of unbound quarks, whose interaction becomes ever weaker with increas¬ 
ing density, approaching the limit of asymptotic freedom predicted by QCD. Any quark can 
now move freely throughout the medium: we have quark liberation through swarm formation. 
Wherever a quark goes, there are many other quarks nearby. The transition from atomic to 
quark matter is schematically illustrated in Fig. [T] 




Figure 1: Schematic view of matter for increasing density, from atomic (a) to nuclear (b) and 
then to quark matter (c). 


We have here considered quark matter formation through the compression of cold nuclear mat¬ 
ter. A similar effect is obtained if we heat a meson gas; with increasing temperature, collisions 
and pair production lead to an ever denser medium of mesons. And according to Hagedom, also 
of ever heavier mesons of an increasing degeneracy. For Hagedom, the fireballs where point¬ 
like, so that the overlap we had just noted simply does not occur. In the real world, however, 
they do have hadronic size, so that they will in fact interpenetrate and overlap before the diver¬ 
gence of the Hagedom resonance gas occurs [fT^ . Hence now again there will be a transition 
from resonance gas to a quark-gluon plasma, now formed by the liberation of the quarks and 
gluons making up the resonances. 

At this point, it seems worthwhile to note an even earlier approach leading to a limiting temper¬ 
ature for hadronic matter. More than a decade before Hagedom, I. Ya. Pomeranchuk IfTTl had 
pointed out that a cmcial feature of hadrons is their size, and hence the density of any hadronic 
medium is limited by volume restriction: each hadron must have its own volume to exist, and 
once the density reaches the dense packing limit, it’s the end for hadronic matter. This simply 
led to a temperature limit, and for an ideal gas of pions of 1 fm radius, the resulting temperature 
was again around 200 MeV. Nevertheless, these early results remained largely unnoticed until 
the work of Hagedom. 

Such geometric considerations do, however, lead even further. If hadrons are allowed to in¬ 
terpenetrate, to overlap, then percolation theory predicts two different states of matter [I181I19II : 
hadronic matter, consisting of isolated hadrons or finite hadronic clusters, and a medium formed 
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as an infinite sized cluster of overlapping hadrons. The transition from one to the other now be¬ 
comes a genuine critical phenomenon, occurring at a critical value of the hadron density. 

We thus conclude that the pioneering work of Rolf Hagedorn opened up the field of critical 
behavior in strong interaction physics, a field in which still today much is determined by his 
ideas. On a more theoretical level, the continuation of such studies was provided by finite 
temperature lattice QCD, and on the more experimental side, by resonance gas analyses of the 
hadron abundances in high energy collisions. In both cases, it was found that the observed 
behavior was essentially that predicted by Hagedorn’s ideas. 


4 Resonance gas and QCD thermodynamics 

With the formulation of Quantum Chromodynamics (QCD) as a theoretical framework for the 
strong interaction force among elementary particles it became clear that the appearance of the 
ultimate Hagedorn temperature Th signals indeed the transition from hadronic phase to a new 
phase of strongly interacting matter, the quark-gluon plasma (QGP) If20l . As QCD is an asymp¬ 
totically free theory, the interaction between quarks and gluons vanishes logarithmically with 
increasing temperature, thus at very high temperatures the QGP effectively behaves like an ideal 
gas of quarks and gluons. 

Today we have detailed information, obtained from numerical calculations in the framework 
of finite temperature lattice Quantum Chromodynamics 11211 - 12511 . about the thermodynamics 
of hot and dense matter. We know the transition temperature to the QGP and the tempera¬ 
ture dependence of basic bulk thermodynamic observables such as the energy density and the 
pressure ll251[26ll . We also begin to have results on fluctuations and correlations of conserved 
charges [lT7] - [2^ . 

The recent increase in numerical accuracy of lattice QCD calculations and their extrapolation 
to the continuum limit make it possible to confront the fundamental results of QCD with Hage- 
dorn’s concepts [|Tl[51, which provide a theoretical scenario for the thermodynamics of strongly 
interacting hadronic matter [|29l - [32l . 

In particular, the equation of state calculated on the lattice at vanishing and finite chemical po¬ 
tential, and restricted to the confined hadronic phase, can be directly compared to that obtained 
from the partition function (fT4l) of the hadron resonance gas, using the form (fT5l) introduced by 
Hagedorn for a continuum mass spectrum. Alternatively, as first approximation, one can also 
consider a discrete mass spectrum which accounts for all experimentally known hadrons and 
resonances. In this case the continuum partition function of the Hagedorn model is expressed 
by Eq. (fT3l) with p(m,) replaced by the spin degeneracy factor of the hadron, and with the 
summation taken over all known resonance species listed by the Particle Data Group If33l . 

With the above assumption on the dynamics and the mass spectrum, the Hagedorn resonance 
gas partition function |[Il[5l can be calculated exactly and expressed as a sum of one-particle 
partition functions Z/ of all hadrons and resonances. 


inz(r,y) = ^z/(r,y). 

i 


(16) 
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For particles of mass m,- and spin degeneracy factor gj, the one-particle partition function Zj, in 
the Boltzmann approximation, reads 


I / X VTmf .nii 


(17) 


Due to the factorization of the partition function in Eq. (fT^ . the energy density and the pressure 
of the Hagedom resonance gas with a discrete mass spectrum, can also be expressed as a sum 
over single particle contributions 


with 


and 


£ = ££/ , p = 


I 

2^4 


Si 


2%'^ ' T 


3i^2(/3m, 

/3m 


■+ 


P} 


Si 


r4 2n^ ' T 


^2 (/3 m,), 


(18) 


(19) 


( 20 ) 


where /3 = l/T and K\ and K 2 are modified Bessel functions. At vanishing chemical potentials 
and at finite temperature, the energy density £, the entropy density s and the pressure P, are 
connected through the thermodynamic relation. 


e^-P + sT. (21) 

Summing up in Eq. (fTSl) the contributions from experimentally known hadronic states, consti¬ 
tutes the resonance gas [fni5ll34l for the thermodynamics of the hadronic phase of QCD. Taking 
into account contributions of all mesonic and baryonic resonances with masses up to 1.8 GeV 
and 2.0 GeV, respectively, amounts to 1026 resonances. 


The crucial question thus is if the equation of state introduced by Hagedorn can describe the 
corresponding results obtained from QCD within lattice approach. In Fig. [2] we show the tem¬ 
perature dependence of the energy density, pressure and the entropy density obtained recently 
in lattice QCD studies with physical masses of up, down and strange quarks [l26l . The bands in 
lattice QCD results indicate error bars due to extrapolation to the continuum limit. The vertical 
band marks the temperature, = (154±9) MeV, which within error, is the crossover temper¬ 
ature from a hadronic phase to a quark-gluon plasma If35ll . These QCD results are compared in 
Fig. |2]to Hagedom’s resonance gas model (HRG) formulated for a discrete mass spectrum in 
Eqs. (fT^ and (|2TI) . There is excellent agreement between the Hagedorn model results for the 
equation of states and the corresponding lattice data. All bulk thermodynamical observables 
are very strongly changing with temperature when approaching the deconfinement transition 
temperature. In Hagedorn’s formulation, this behavior is well understood in terms of increasing 
resonance contributions. Although the HRG formulated for discrete mass spectrum does not 
exhibit any critical behavior, it nevertheless reproduces remarkably well the lattice results in 
the hadronic phase. This agreement has now been extended to an analysis of fluctuations and 
correlations of conserved charges as well. 
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Figure 2: The normalized pressure /’(T), the energy density e(r) and the entropy density s{T) 
obtained in lattiee QCD ealeulations as a funetion of temperature. The dark lines show the 
predietion of the Hagedorn resonanee gas for diserete mass speetrum, eqs. (fT^ - (l2TI) . The 
lattiee results are from Ref. If26ll . 

The exeellent deseription of the lattiee QCD results by Hagedorn’s model justifies the elaim that 
resonances are indeed the essential degrees of freedom near deconfinement. Thus, on the ther- 
modynamieal level, modeling hadronie interaetions by formation and exeitation of resonanees, 
as introdueed by Hagedorn, is an exeellent approximation of strong interaetions. 


5 Resonance Gas and Heavy Ion Collisions 


Long before lattiee QCD eould provide a direet evidenee that strong interaetion thermodynam- 
ies ean be quantified by the resonanee gas partition funetion, Hagedorn’s eoneept was veri¬ 
fied phenomenologieally by eonsidering partiele produetion in elementary and heavy ion eolli- 
sions Il^l36l - [39ll . In a strongly interaeting medium, one ineludes the eonservation of eleetrie 
eharge, baryon number and strangeness. In this ease, the partition funetion of Hagedom’s ther¬ 
mal model depends not only on temperature but also on ehemieal potential /I, whieh guarantees, 
that eharges are eonserved on an average. For a non vanishing /I, the partition funetion Eq. (fT^ 
is replaeed by 

inZ(r,y,M) = i;z/(r,y,At), (22) 

i 

with = (hb-iIJ-s-iIJ-q), where /i,- are the ehemieal potentials related to the baryon number, 
strangeness and eleetrie eharge eonservation, respeetively. 

For partiele i earrying strangeness Si, the baryon number Bi, the eleetrie eharge Qi and the 
spin-isospin degeneraey faetor gt, the one partiele partition funetion, reads 


zi{T,v,n) 


VgiTmj 


^2(w,/r)exp 


^il^B -fSilAs + QiBQ ^ 


(23) 


For p. = 0 one reeovers the result from Eq. (fTV]) . 
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The calculation of a density ni of particle i from the partition function Eq. (I22l l is rather straight¬ 
forward [[soil . It amounts to the replacement Z/ —)■ jiZ] in Eq. (l22l) and taking a derivative with 
respect to the particle fugacity ji, as 


rii = 


m 


th 


V 


dlnZ 


dji 


Yi=l 


(24) 


consequently, n, = Z/ jW with Z\ as in Eq. (|2^ . 

The Hagedom model, formulated in Eq. (l22l) . describes bulk thermodynamic properties and 
particle composition of a thermal fireball at finite temperature and at non vanishing charge 
densities. If such a fireball is created in high energy heavy ion collisions, then yields of different 
hadron species are fully quantified by thermal parameters. However, following Hagedorn’s idea, 
the contribution of resonances decaying into lighter particles, must be included [UllSlI. 

In Hagedom’s thermal model, the average number (A^,) of particles i in volume V and at temper¬ 
ature r that carries strangeness 5/, the baryon number Bu and the electric charge Qu is obtained 
from Eq. (l22l) . see [[BEl 


(w.Xr.fi) = + (25) 

The first term in Eq. (|251) describes the thermal average number of particles of species i from Eq. 
(I24l) and the second term describes overall contribution from resonances. This term is taken as a 
sum of all resonances that decay into particle i. The T is the corresponding decay branching 
ratio of j —?■ i. The multiplicities of resonances {NjY^’^ in Eq. (1251) . are obtained from Eq. (I24l) . 

The importance of resonance contributions to the total particle yield in Eq. (1251) is illustrated 
in Eig. dU) for charge pions. In Eig. ([3]) we show the ratio of the total number of charge pions 
from Eq. (1251) and the number of prompt pions from Eq. (l24l) . The ratio is strongly increasing 
with temperature and chemical potential. This is due to an increasing contribution of mesonic 
and baryonic resonances. Erom Eq. ([3]) it is clear, that at high temperature and/or density, the 
overall multiplicity of pions is mostly due to resonance decays. 

The particle yields in Hagedorn’s model Eq. (l25l) depend, in general, on five parameters. How¬ 
ever, in high energy heavy ion collisions, only three parameters are independent. The isospin 
asymmetry, in the initial state fixes the charge chemical potential and the strangeness neutrality 
condition eliminates the strange chemical potential. Thus, on the level of particle multiplicity, 
we are left with temperature T and the baryon chemical potential jig as independent parameters, 
as well as, with fireball volume as an overall normalization factor. 

Hagedorn’s thermal model introduced in Eq. (l25l) was successfully applied to describe particle 
yields measured in heavy ion collisions. The model was compared with available experimental 
data obtained in a broad energy range from AGS up to EHC. Hadron multiplicities ranging 
from pions to omega baryons and their ratios, as well as composite objects like e.g. deuteron or 
alpha particle, were used to verify if there was a set of thermal parameters (T, /ig) and V, which 
simultaneously reproduces all measured yields. 

The systematic studies of particle production extended over more than two decades, using ex¬ 
perimental results at different beam energies, have revealed a clear justification, that in central 
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Figure 3: The ratio of the total density of positively eharged pions, from Eq. (l25l) . and 
the density of thermal pions, from Eq. (l24l) . The ealeulations are done in the Hagedom 
resonanee gas model for =250 MeV and Hb =550 MeV at different temperatures. 


heavy ion eollisions partiele yields are indeed eonsistent with the expeetation of the Hagedom 
thermal model. There is also a elear pattern of the energy, i/y-dependenee of thermal param¬ 
eters. The temperature is inereasing with ^/s, and at the SPS energy essentially saturates at 
the value, whieh eorresponds to the transition temperature from a hadronie phase to a QGP, as 
obtained in EQCD. The ehemieal potential, on the other hand, is gradually deereasing with ^/s 
and almost vanishes at the EHC. 


In Pig. dH) we show, as an illustration, a eomparison of Hagedom’s thermal model and reeent 
data on seleeted partiele yields, obtained by AEICE eollaboration in eentral Pb-Pb eollisions 
at midrapidity at the EHC energy [l40ll4Tll . At sueh high eollision energy, partiele yields from 
Eq. (1251) are quantified entirely by the temperature and the fireball volumeQ Thus, there is 
transparent predietion of Hagedom’s model Eq. (l25l) . that yields of heavier partieles {Ni) with 
no resonanee deeay eontributions, normalized to their spin degeneraey faetor gi = (27 -f 1), 
should be quantified by 


W-) 

Z/+ 1 


Q / nii \3/2 

\2nT 


j exp(-m//r), 


(26) 


where we have used Eq. (l24l) and the asymptotie expansion of the Bessel funetion, 7^2 (x) ~ 
ji;^^/^exp(—x), valid for large x. 

In Pig. dH) we show the yields of partieles with no resonanee eontribution, like (j), Q, the 
deuteron ‘d’, ^He and the hypertriton ^He, normalized to their spin degeneraey faetor, as a 
funetion of partiele mass. Also shown in this figure is the predietion from Eq. (|2^ at T ~ 156 
and for volume V ~ 5000 fm^ HOlldTll . There is a elear eoineidenee of data taken in Pb-Pb 
eollisions at the EHC and predietions of the Hagedom model Eq. (|2^ . Partieles with no 
resonanee eontribution measured by AEICE eollaboration follow the Hagedom’s expeetations 
that they are produeed from a thermal fireball at eommon temperature. A similar agreement 

*The chemical potential jl in Eq. (l25l) vanishes, since at the LHC and at midrapidity particles and their antipar¬ 
ticles are produced symmetrically. 
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Figure 4: Yields of several different particle species per unit rapidity normalized to spin de¬ 
generacy factor as a function of their mass. Data are from ALICE collaboration taken at the 
LHC in central Pb-Pb collisions. The line is the Hagedom thermal model result, Eq. (|2^ . see 
Ref. [|40ll4T1. 

of Hagedom’s thermal concept and experimental data taken in central heavy ion collisions has 
been found for different yields of measured particles and collision energies from AGS, SPS, 
RHIC and LHC [|3l- 

6 Particle Yields and Canonical Charge Conservation 

The Hagedom thermodynamical model for particle production, was originally applied to quan¬ 
tify and understand particle yields and spectra measured in elementary collisions - there were 
no data available from heavy ion collisions. 

Initial work on particle production by Hagedom begins in 1957 in collaboration with E. Cerulus 
when they apply the Eermi phase space model. In this microcanonical approach, conservation 
laws of baryon number or electric charge are implemented exactly. Almost 15 years later the 
production of complex light antinuclei, such as anti-He^, preoccupied Hagedom [[IllSl. He re¬ 
alized and discussed clearly the need to find a path to enforce exact conservation of baryon 
number to describe the anti-He^ production correctly within the canonical statistical formula¬ 
tion. 

Indeed, applying in pp reactions the thermal model without concern for conservation of baryon 
number overestimates the production of anti-He^ in proton-proton collisions by seven orders 
of magnitude [|Tl|3l[5l. The reason was that when the number of particles in the interaction 
volume is small, one has to take into account the fact that the production of anti-He^ must be 
accompanied by the production of another three nucleons with energy E;^, in order to exactly 
conserve the baryon number. Thus, in case that the production of anti-He^ is not originating 
from reservoir of many antiquarks or antinucleons already present in a large volume, but is 
rather originating from some small volume Vpp that is present in pp collisions, the abundance 
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of anti-He^ will not be proportional to the single standard Boltzmann factor, as in Eq. (|2^ 


~ exp 



( 27 ) 


but is accompanied by additional Boltzmann factors that characterize the production of the 
associated nucleons, needed in order to conserve baryon number [dllSl 



This suppresses the rate and introduces a strong power-law dependence on volume Vpp for the 


anti-He^ yield. 


The problem of exact conservation of discrete quantum numbers in a thermal model formulated 
in early ’70s by Hagedom in the context of baryon number conservation remained unsolved for a 
decade. When the heavy ion QGP research program was approaching and strangeness emerged 
as a potential QGP signature, Hagedorn pointed out the need to consider exact conservation 
of strangeness ll^ . This is the reason that the old problem of baryon number conservation 
was solved in the new context of strangeness conservation [|43T - |45]| . A more general solution, 
applicable to all discrete conserved charges, abelian and non-abelian, was also introduced in 
Ref. ||46]| and expanded in [|T7] - l52l . Recently, it has become clear that a similar treatment 
should be followed not only for strangeness but also for charm abundance study in high energy 
collisions II531I54II . 

To summarize this section, we note that the usual form of the statistical model, based on a grand 
canonical formulation of the conservation laws, cannot be used when either the temperature or 
the volume or both are small. As a rough estimate, one needs VT^ > I for a grand canonical de¬ 
scription to hold [I431I50II . In the opposite limit, a path was found within the canonical ensemble 
to enforce charge conservations exactly. 

The canonical approach has been shown to provide a consistent description of particle produc¬ 
tion in high energy hadron-hadron, and peripheral heavy ion collisions |IMl|49l[53l[54l . In 
the context of developing strangeness as signature of QGP, such a model also provides, within 
the realm of assumed strangeness chemical equilibrium, a description of an observed increase 
of single- and multi-strange particle yields from pp, pA to AA collisions and its energy depen¬ 
dence Il44l . 

7 Concluding Remarks 

Rolf Hagedorn’s work, introducing concepts from statistical mechanics and from the mathe¬ 
matics of self-similarity into the analysis of high energy multiparticle production, started a new 
field of research, alive and active until today. On the theory side, the limiting temperature of 
hadronic matter and the behavior of the Hagedorn resonance gas approaching that limit were 
subsequently verified by first principle calculations in finite temperature QCD. On the experi¬ 
mental side, particle yields as well as, more recently, fluctuations of conserved quantities, were 
also found to follow the pattern predicted by the Hagedom resonance gas. Rarely has an idea in 
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physics risen from such humble and little appreeiated beginnings to sueh a striking vindieation. 
So perhaps it is appropriate to close with a poetie summary one of us (HS) formulated some 
twenty years ago for a Hagedorn-Fest, with a slight update. 


Hot Hadronic Matter 

(A Poetic Summary) 

In days of old 
a tale was told 
of hadrons ever fatter. 
Behold, my friends, said Hagedorn, 
the ultimate of matter. 

Then Muster Mark 
called in the quarks, 
to hadrons they were mated. 

Of eolors three, and never free, 
all to confinement fated. 

But in dense matter, 
their bonds ean shatter 
and they can freely move around. 
Above Tfi, their colors shine, 
as the QGP is found. 

Said Hagedorn, 
when quarks were born 
they had different advances. 
Today they form, as we ean see, 
a gas of all their chances. 
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